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THE EFFECTS OF DUST IN SIMPLE ENVIRONMENTS: LARGE MAGELLANIC CLOUD H II REGIONS 
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ABSTRACT 

We investigate the effects of dust on Large Magellanic Cloud (LMC) H II region spectral energy distributions 
using arcminute-resolution far-ultraviolet (FUV), Ha, far-infrared (FIR), and radio images. Widely-used indica- 
tors of the amount of light lost to dust (attenuation) at Ha and in the FUV correlate with each other, although often 
with substantial scatter. There are two interesting systematic discrepancies. First, Ha attenuations estimated from 
the Balmer decrement are lower than those estimated from the Ha-to-thermal radio luminosity ratio. Our data, 
at this stage, cannot unambiguously identify the source of this discrepancy. Second, the attenuation at 1500A 
and UV spectral slope, j3, correlate, although the slope and scatter are substantially different from the correla- 
tion first derived for starbursting galaxies by Calzetti et al. Combining our result with those of Meurer et al. for 
ultra-luminous infrared galaxies and Calzetti et al. for starbursting galaxies, we conclude that no single relation 
between (3 and 1500A attenuation is applicable to all star-forming systems. 

Subject headings: dust, extinction — H II regions — Magellanic Clouds — galaxies: individual (LMC) — 
ultraviolet: galaxies 

1. INTRODUCTION 



Estimating where and how quickly stars are forming in galax- 
ies of different types, as a function of redshift, gives us direct in- 
sight into where, when and how galaxies evolve (see e.g. Madau 
et al. 1996; Blain et al. 1999; Steidel et al. 1999, for attempts to 
address this problem in a global sense). In this context, under- 
standing the systematic uncertainties plaguing SFR estimators 
is of special importance. One of the principal uncertainties of 
this type is the correction for emergent flux lost due to dust. In a 
systematic effort to address this problem, we begin by studying 
the effects of dust in the simplest of all star formation environ- 
ments, H II regions. 

Before discussing the effects of dust, it is necessary to un- 
derstand the nature and origins of the different SFR indicators. 
SFRs are usually estimated from spectral regions where the en- 
ergy output from young stars dominates. Massive, short-lived 
stars produce copious quantities of ultraviolet (UV) radiation. 
This UV light may be directly observed, or dust may absorb 
much of it, reprocessing it into the far-infrared (FIR). In the 
limit of dust absorbing all of the energy output of young stellar 
populations, the FIR radiation is approximated by the bolomet- 
ric luminosity. The time evolution of the UV and bolometric lu- 
minosities is shown in Fig. 1 for a single burst from the PEGASE 
stellar population models (Fioc & Rocca-Volmerange 2001): 
the UV and bolometric luminosities evolve similarly for more 
than 100 Myr. The most massive and short-lived stars produce 
radiation capable of ionizing the hydrogen gas in their natal 
clouds, producing a H II region. This ionized hydrogen plasma 
produces line emission (the Balmer series is the most frequently 
observed series of emission lines) and thermal bremsstrahlung 
in the radio: these radiations evolve much more quickly than 
the 1500A UV or bolometric luminosities (Fig. 1). 

However, the Ha and UV luminosities will be strongly re- 
duced by dust. The amount of light lost from a system is re- 
ferred to as the attenuation. The attenuation and the more com- 
monly used extinction are the same for a single star with dust 
along the line-of-sight. However, for more complex systems, 
the attenuation is a function of the geometry of both the dust 
and stars. Classic methods that work for determining extinc- 
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FIG. 1. — Age and metallicity dependence of the bolometric luminosity 
(Lboi ; solid lines), 1500A luminosity (L1500; dotted lines), and Ha luminos- 
ity (Lhq ; derived assuming 0.45 Ha photons per Lyman continuum photon; 
dashed lines) of a IMq stellar population, as derived using the PEGASE stel- 
lar population synthesis models (Fioc & Rocca-Volmerange 2001). Thick lines 
denote the Z=0.008 (40% solar) models, thin lines denote the Z=0.004 (20% 
solar) models, and thin lines with crosses denote the Z=0.02 (solar) models. 
All luminosities at ages > 3 Myr decrease with increasing metallicity, primar- 
ily due to increased stellar opacities. At younger ages, the behavior is more 
complex because of the interplay of mean molecular weight and opacity ef- 
fects. 



tions (e.g. color excesses) will not reliably determine attenua- 
tions. 

Despite the difficulties introduced by geometry, it is possible 
to estimate the amount of luminosity lost at Ha or in the UV, al- 
lowing the true SFR to be estimated. Some attenuation estima- 
tors rely solely on UV/optical data. For example, it is common 
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to estimate the amount of Ha light lost to dust by comparing 
the Ha and H/3 emission-line fluxes, which have an intrinsic 
ratio fixed by atomic physics in most astrophysical environ- 
ments. The observed relative ratio therefore gives an estimate 
of the reddening: this can be transformed into an attenuation 
under the assumption of a particular attenuation curve (see e.g. 
Quillen & Yukita 2001, for a test of this method, and references 
therein). A common approach to estimate the attenuation of UV 
stellar light is to use the UV spectral slope (3. This approach, 
pioneered by Calzetti et al. (1994) and Meurer, Heckman, & 
Calzetti (1999), uses the relatively narrow range in intrinsic (3 
for young, luminous stellar populations (-2.6 < (3 < -1.5 in- 
cluding even metal-rich single bursts with ages younger than 
20 Myr; Leitherer et al. 1999) to estimate the reddening. This 
reddening is observed to correlate with UV attenuation with 
only modest scatter for local starbursting galaxies (Calzetti et 
al. 1994; Heckman et al. 1998; Meurer, Heckman, & Calzetti 
1999). However, if one applies this relationship to systems 
other than those for which it was calibrated, then one is as- 
suming that the relation holds over a variety of star-forming 
systems. 

Other attenuation estimators rely on data from beyond UV 
and optical wavelengths. For example, the Ha attenuation can 
be estimated from the ratio of Ha to thermal radio contin- 
uum (e.g. Israel & Kennicutt 1980; Kennicutt 1983; van der 
Hulst et al. 1988). In this case, one concern is that the ther- 
mal radio emission may be contaminated by non-thermal emis- 
sion. This contamination is severe for most galaxies (in fact, 
the non-thermal radio may correlate with SFR, although we do 
not discuss it in this paper; see e.g. Condon 1992; Cram et al. 
1998). Another such attenuation indicator, this time for the far- 
UV (FUV), uses the ratio of FUV to FIR, which measures the 
amount of FUV light reprocessed into the FIR (e.g. Buat & Xu 
1996; Meurer, Heckman, & Calzetti 1999; Gordon et al. 2000). 

There are inevitable limitations to any attenuation estimator. 
Bearing in mind the importance of constructing relatively ac- 
curate, unbiased SFRs as a probe of galactic evolution, it is 
vital to critically assess the applicability and efficacy of atten- 
uation estimators. Considering data from as many wavelengths 
as possible is an important part of this assessment: by using 
UV, Ha, FIR, and thermal radio luminosities it is possible to 
construct different estimates of the attenuation at a given wave- 
length. Examples of this kind of multi-wavelength approach 
for galaxies include Calzetti et al. (1994), Meurer, Heckman, 
& Calzetti (1999), Sullivan et al. (2000), and Bell & Kennicutt 
(2001). 

An essential test of the robustness of attenuation indicators 
is their behavior in the simplest of all composite star-forming 
systems: Hll regions. Without a firm understanding of the 
attenuation of Ha and FUV light in H II regions, we cannot 
claim a firm understanding of the effects of dust on galactic 
scales. Because Hll regions have ionized gas and relatively 
well-constrained ages (the dominant populations have to be ca- 
pable of ionizing hydrogen), it is possible to intercompare the 
attenuations derived for the ionized gas and those derived for 
the stellar population inside the H II region. We have chosen to 
use the sample of Large Magellanic Cloud (LMC) H II regions 
from Caplan & Deharveng (1986) for this purpose. This sam- 
ple of H II regions has a number of advantages for this purpose: 
a common distance, modest galactic foreground extinction, and 
the availability of large-format images and published data in the 
UV, optical, FIR, and radio spectral regions with a useful physi- 



cal resolution (we adopt 4.9' [80 pc] diameter apertures for this 
work). 

The remainder of this paper is organized as follows. §2 dis- 
cusses the construction and limitations of FUV and Ha atten- 
uation indicators. §3 describes the multi-wavelength imaging 
data, the morphology of the LMC at different wavelengths, and 
the photometry of individual Hll regions. §4 quantifies the 
relationships between different H II region attenuation indica- 
tors. In §5, we discuss in detail two particular comparisons: 
the well-documented discrepancy between Ha attenuations es- 
timated from Balmer decrements and Ha-to-radio ratios, and 
the highly scattered and rather steep correlation between FUV 
attenuation and UV spectral slope (3. §6 summarizes our con- 
clusions from this study. We assume a distance to the LMC 
of 45.9 ±2 kpc (e.g. Fitzpatrick et al. 2001), which only af- 
fects the absolute scaling of the SEDs presented in §3 because 
the attenuations and flux ratios are distance-independent. In 
correcting our attenuation estimates for foreground galactic ex- 
tinction, we adopt a foreground reddening of E(B-V) = 0.06, 
following Schlegel et al. (1998) and Oestreicher, Gochermann, 
&Schmidt-Kaler(1995). 

2. QUANTITATIVE ATTENUATION INDICATORS 

In this section, we discuss the estimation of Ha and UV at- 
tenuations from multiwavelength data. This serves a dual pur- 
pose: to introduce the methods used to calculate the attenuation 
indicators later in the paper, and to instill an intuition for what 
to expect when comparing the different attenuation indicators, 
including the nontrivial effects of radiative transfer. 

2.1. Attenuation of far- UV light 

Using measured UV reddenings to estimate the FUV attenu- 
ation is becoming increasingly fashionable due to its observa- 
tional efficiency. The UV continuum spectral slope (3, defined 
by: 

Fx cx \P, (1) 

where F\ is the flux per unit wavelength A, is used to estimate 
the attenuation in the UV (e.g. Calzetti et al. 1994; Meurer, 
Heckman, & Calzetti 1999; Adelberger & Steidel 2000) be- 
cause one presumes to understand the color of the underlying 
starlight, and so any deviation from this color is attributed to 
dust. To the degree that one believes the underlying assump- 
tions in this method, one could, for example, also use optical 
reddenings to estimate UV attenuation, or use UV reddenings 
to estimate the optical attenuation. Clearly, such an exagger- 
ated extension of the technique loses credibility and one begins 
to wonder where, exactly, the technique may break down. The 
relationship between the UV or optical reddenings and the FUV 
attenuation is non-trivial, as illustrated in panels (a) and (b) of 
Fig. 2. 

In Fig. 2, we show predictions from the radiative transfer 
model DIRTY (Witt & Gordon 2000; Gordon et al. 2001; Mis- 
selt et al. 2001). The relationships between different attenuation 
indicators are shown under different assumptions about the ex- 
tinction curve (MW=Milky Way, LMC2=a supershell near 30 
Dor, SMC=Small Magellanic Cloud Bar), dust clumping (ho- 
mogeneous or clumpy, with a density ratio of 100), and geom- 
etry (dusty=uniformly mixed stars and dust, shell=inner dust- 
free region of stars surrounded by a shell of star- free dust). The 
main point to note from the panels (a) and (b) of Fig. 2 is that it 
is possible to have a large range in FUV attenuation for a given 
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FIG. 2. — Predictions from radiative transfer models regarding the relationships between different attenuation indicators, (a) 1500A attenuation against change in 
UV spectral slope A/3, (b) E(B-V) against 1500A attenuation, (c) Balmer-derived Ha attenuation against true Ha attenuation, (d) B aimer-derived Ha attenuation 
against 1500A attenuation, and (e) true Ha attenuation against 1500 A attenuation. Models are described in the legend, and expanded on in the text. Panel (a) also 
shows screen models with LMC (solid line) and SMC Bar (dashed line) extinction curves. 



UV or optical reddening, depending on the dust/star geometry 
and/or shape of extinction curve. Bearing in mind the exten- 
sive empirical evidence for a highly environmentally-sensitive 
UV extinction curve (e.g. Fitzpatrick 1986; Gordon & Clayton 
1998; Misselt, Clayton & Gordon 1999; Clayton, Gordon, & 
Wolff 2000), and the diversity of dust/star geometries observed 
in many astrophysical systems, our expectation is that the va- 
lidity of UV or optical reddening-based attenuations is open to 
suspicion in any but the best-studied and characterized cases. 

It is possible to avoid the need to assume a particular extinc- 
tion curve in deriving the FUV attenuation by comparing the 
observed FUV starlight with the FIR luminosity. Because the 
majority of the bolometric flux of a young stellar population 
is emitted in the FUV and because the FIR flux is dominated 
by the re-radiation of attenuated UV light, this comparison pro- 
vides an estimate of the FUV attenuation that is insensitive to 
the extinction curve and geometry when averaged over 4ir stera- 



dians (see Figs. 1 and 2 of Gordon et al. 2000). We use an 
IDL subroutine, available from K. D. G, and described fully in 
Gordon et al. (2000), to estimate the attenuation at 1500A and 
1900A from the 1500A, 1900A, and FIR fluxes. This subrou- 
tine takes into account the full uncertainties implied by differ- 
ences in dust geometry, dust type and stellar populations. For 
1500 A attenuations between and 4 mags, the 1500 A flux ra- 
tio method-derived attenuation can be approximated by a mod- 
ification of the simple energy-balance estimate to within 0.03 
mag: 

Ai5oo,fr~0.86 -2.51og 10 — ^ 0.1, (2) 

VJv + i-FIR _ 

where the bracketed quantity is the simple energy-balance esti- 
mate of 1500A attenuation {vf v is evaluated at 1500A), Lfir is 
an estimate of the total FIR luminosity (estimated as described 
in §4.2.1), and the scaling and subtraction approximately cor- 
rect for SED and attenuation curve shape. This method has 
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important limitations. First, it requires isotropy on scales larger 
than the illuminating source. An example of a geometry for 
which this approach fails is the case of a torus of dust outside a 
cluster of FUV bright stars, in which the torus is oriented per- 
pendicular to the line of sight. Such an region would have a 
large FIR luminosity but no attenuation along the line of sight. 
Second, the method requires the young population to dominate 
the heating. The method, some of its uncertainties, and its lim- 
itations are presented in Gordon et al. (2000): variants of this 
method were used by e.g. Buat & Xu (1996) or Meurer, Heck- 
man, & Calzetti (1999) to explore the attenuations of spiral and 
starburst galaxies respectively. 

2.2. Attenuation of Ha light 

Estimating the Ha luminosity lost to dust is relatively straight- 
forward. We discuss two methods: both estimate only the at- 
tenuation of the Ha radiation and are insensitive to the amount 
of Lyman continuum radiation lost to dust before being down- 
converted into Ha. The first is based on the comparison of the 
Ha line luminosity to the thermal radio luminosity. This ratio, 
in the absence of dust, is given by equations A5 and A10 of 
Caplan & Deharveng (1986): 

/h q //8.55GHz = 8.659 x 10" 9 ergcm" 2 s _1 Jy" 1 

x (t e )-° M [8 .666 + 1 .5 lnfc)] -1 , (3) 

where fu a is the Ha flux in erg cm" 2 s" 1 , /s.ssghz is the 8.55 
GHz flux in Jy, and t e is the electron temperature in units of 
10 4 K. The ratio depends weakly on electron temperature: we 
will adopt an electron temperature of 9500 K ± 1000K (Pagel 
et al. 1978), corresponding in dust-free Ha to radio ratio un- 
certainties of ~8%. In practice, the thermal radio luminosity 
is used to construct an expectation for the Ha line luminos- 
ity. Then, the observed Ha line luminosity is compared to the 
radio expectation to derive an effective attenuation Ana.Radio- 
Although this should be a robust attenuation indicator, accurate 
thermal radio luminosities are difficult to measure, even for in- 
dividual H II regions. 

The second approach is to estimate the Ha attenuation using 
the observed Ha/H/3 ratio for a given H II region. Assuming a 
standard Ry = 3 . 1 Milky Way extinction curve and a foreground 
screen dust distribution, Caplan & Deharveng (1986) derive 
the following expression for Ha attenuation estimated from the 
Balmer decrement A Ha ,B a imer given an observed Ha/H/3: 

/ Ha/H/3 \ 

A HQ ,Balmer = 5.251og 10 I - ggg^xn ) ' ( 4 ) 

or, equivalently, A Ha B . dlmer = 2.1 E(Hf3- Ha). Ha/H/3 is obser- 
vationally straightforward to measure; however, different dust 
extinction curves or geometries change the conversion between 
reddening and attenuation from the standard R v = 3.1 fore- 
ground screen. 

The sensitivity of these estimators to geometry and dust prop- 
erties is illustrated in panel (c) of Fig. 2, where A H a,B a imer is 
shown against the true Ha attenuation A(Ha). For homoge- 
neous shell or screen geometries, AHa.Baimer and A(Ha) are es- 
sentially identical; however, clumpy shell and dusty geometries 
show a 'saturation' for A Ha ,Baimer ^0.8 and ^0.4 mag respec- 
tively, where only the least reddened stars are observed, result- 
ing in high attenuations but modest reddenings. 

2.3. Comparing Ha andfar-UV attenuations 



The relationship between Ha and FUV attenuations is com- 
plex because it depends on extinction curve shape, dust clump- 
ing properties, and geometry (panel (d) of Fig. 2). In particular, 
for clumped dust or mixed star/dust geometries, the relationship 
between AHa.Baimer and 1500A attenuation is non-linear. Panel 
(d) is quite similar to panel (b); both A HQ ,B a imer and E(B-V) 
are essentially color excesses, and behave quite similarly. Also, 
as shown in panel (e), the relationship between true Ha and 
1500 A attenuation is non-linear for clumped and/or dusty ge- 
ometries (in quite a different way than the Balmer-derived Ha 
attenuation). While the FUV attenuation 'saturates' at around 
2 magnitudes, the Ha attenuation continues to increase (be- 
cause the FUV becomes optically thick, but the Ha remains 
optically thin). A further complication, which is more of a con- 
cern when comparing Ha and FUV attenuation in galaxies, is 
that the Ha and FUV light can come from spatially distinct lo- 
cations. In H II regions, this is less of a concern, and we will 
use the inter-comparison of Ha and FUV attenuations to test 
the reproducibility of these attenuation estimators. 

3. DATA ANALYSIS 

3.1. Images 

We use FUV, Ha, FIR, and 8.55 GHz radio images of the 
LMC to investigate the spectral energy distributions (SEDs) and 
attenuations for H II regions drawn from the catalogs of Caplan 
& Deharveng (1985, 1986). Accordingly, we have compiled 
an extensive library of LMC images with roughly arcminute 
resolution from the literature. We describe these images in turn 
below, and in Table 1 . 

• The UV images come from a rocket launched from Aus- 
tralia in 1977 (for a full description see Smith, Cornett, & Hill 
(1987)). The angular resolution is ^100" over most of the use- 
ful part of the image. We convert the raw digital units to MJy/sr 
using calibrations supplied by R. Cornett (1999, private com- 
munication). We rule out offsets in calibration larger than 12% 
at the 2a level by comparing to aperture fluxes of 10 star form- 



Table 1 
LMC Data 



Band 


A 


AA FWHM 


Ref 




[/xm] 


[H 


["] 




RocketJJVl 


0.1495 


0.020 


100 


1 


Rocket_UV2 


0.1934 


0.022 


100 


1 


R 


0.66 


0.15 


72 


2 


Ha 


0.6571 


0.0013 


72 


2 


Msxm 


8.3 


3.36 


45 


3 


IRAS 12,um 


12 


6.5 


38 


4 


IRAS 25,um 


25 


11 


42 


4 


IRAS 60/zm 


60 


40 


75 


4 


IRAS lOO^rn 


100 


37 


135 


4 


8550 MHz 


3.5 x 10 4 




210 


5 
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Hill (1987), (2) Kennicutt et al. (1995), (3) 
http://www.ipac.caltech.edu/ipac/msx/msx.html, (4) Braun 
et al. (1998), (5) Filipovic et al. (1998) 
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ing regions from Ultraviolet Imaging Telescope images with an 
effective wavelength of 1614A. 

• The Ha and R band images were taken with the Parking Lot 
Camera (Bothun & Thompson 1986), but with a Texas Instru- 
ments 800x800 pixel CCD (Kennicutt et al. 1995). We cali- 
brated the Ha image by comparing to aperture Ha fluxes from 
Caplan & Deharveng (1985). Our calibration is consistent with 
the calibration of the Ha image given by Kennicutt et al. (1995) 
to within 10%. The R band image has not been photometrically 
calibrated. 

• The Midcourse Space Experiment (MSX) B3 (8.3/mi) data 
were downloaded from the IPAC MSX site (for a description 
of the instrument and the data see Price et al. 2001). We use 
the mosaic image, with 45" instead of 18" resolution, as our 
work is limited by the lower resolutions of the images in other 
passbands. 

• The IRAS images at 12/im, 25/«n, 60/im and 100/im, pro- 
cessed with a maximum entropy algorithm to improve resolu- 
tion (Braun et al. 1998), were provided by M. Braun (2000, 
private communication). We confirm the validity of structures 
seen in the resolution-enhanced images by comparing the pro- 
cessed IRAS 12/im image with the MSX 8.3/im image. The 
morphologies of the structures in the images are very similar. 
At faint levels, the processed image displays the characteristic 
striping of IRAS data. 

• The radio images were provided by M. Filipovic (2000, pri- 
vate communication) and are described in detail in Filipovic et 
al. (1998). These data are calibrated to better than 15%, when 
tested against 'aperture' 5 GHz radio fluxes from Caplan & 
Deharveng (1986): before comparison, the 5 GHz radio fluxes 
were scaled down by 5% to account for wavelength differences, 
assuming purely thermal radio emission. 

All the images were aligned (translation and rotation only) to 
the R band image, for which J. Parker (1999, private commu- 
nication) provided an astrometric solution. At least 5-10 stars 
and H II regions were used as the basis of the alignment, and 
alignment was accurate to better than 20" in all cases. After 
the alignment, the images were transformed to a common NCP 
projection with 30" pixels. All of the images, except for the 
Ha and R band image, are calibrated in units of MJy/sr. The 
Ha image is calibrated in ergs s" 1 cm" 2 and, as stated before, 
the R band image is uncalibrated. The original and calibrated, 
aligned images are available from the authors. 

3.2. Morphology of the LMC at different wavelengths 

Before discussing the SEDs and attenuations of the LMC H II 
regions, we must comment on one of the most striking aspects 
of Fig. 3 — the strong difference between the distributions of 
the young stars (as probed by the Ha, 1900 A or FIR emission) 
and older stars (as probed by the R band). The older stellar 
populations are strongly concentrated in the bar, with an under- 
lying disk-like envelope. In contrast, the younger stars have a 
clumpy distribution, and are predominantly concentrated in the 
northern half of the LMC. This interesting aspect of the data is 
beyond the scope of the current work. 

More subtle, but perhaps equally significant, differences are 
also present. In §1, we saw that the UV and FIR luminosities 
should track each other quite closely if the amount and distribu- 
tion of dust around young stars does not depend on stellar age. 
Therefore, a naive expectation might be that the morphologies 
of the LMC in the UV and FIR should be similar. In contrast, 
the Ha morphology could either be similar or not depending 



on the detailed star formation history over timescales of ^ 10 to 
100 Myr. Surprisingly, the observation is that the FIR and Ha 
morphologies are similar (cf. Fig. 3), and that they differ from 
the morphology in the UV. 

This departure from the naive model has two interesting im- 
plications. Firstly, the assumption used in linking UV to FIR 
luminosities has broken down. The similarity between the Ha 
and FIR morphologies implies that there is a good correspon- 
dence between the distribution of very young (ionizing) stars 
and dust. However, as stars age (so that they are now no longer 
strongly ionizing stars but still have relatively large UV fluxes), 
the distribution of dust must decouple from that of the stars. 
We suggest the simple scenario in which young, dusty H II re- 
gions are dispersed by the cumulative effects of stellar winds 
and supernovae. This dispersal (1) causes many of the FUV- 
bright regions to be faint in the FIR, (2) explains the existence 
of large, extended FUV -bright regions ringed by Ha and FIR 
emission from recent star formation, that are triggered by the 
winds from the OB association in the center (e.g. Constellation 
III, in the top left quadrant of the images), and (3) is consistent 
with other empirical evidence (e.g. Chariot & Fall 2000, and 
references therein). 

Secondly, the dissimilarity between the Ha and UV mor- 
phologies of the LMC stresses a fundamental limitation of this 
work in its application to the global situation in galaxies. Be- 
cause much of the UV flux of the LMC comes from outside the 
H II regions we study here, we cannot meaningfully constrain 
the effects of dust on the total UV emission from the LMC (we 
will address this issue in future work). The dichotomy between 
the sites of Ha and UV emission also leads us to caution against 
attempts to link the effects of dust on galactic Ha and UV emis- 
sion too closely: it may be impossible to meaningfully link the 
Ha and FUV extinction properties of star-forming galaxies in 
any more than a statistical sense. 

3.3. SEDs 

We selected a sample of H II regions with measured Balmer 
decrements in 4.9' diameter apertures from the sample of Ca- 
plan & Deharveng (1985, 1986). We retained the use of the 4.9' 
aperture to ensure a common aperture for all attenuations: the 
adoption of a larger aperture (up to 20', depending on source 
size and morphology) gave very similar attenuations at better 
than the 10% level (and similar values of UV spectral slope (3 at 
better than the 0.5 level). The locations of the H II regions are 
indicated in Fig. 3. An explicit selection was made for isola- 
tion: any Hll regions from Caplan & Deharveng (1985, 1986) 
with significant contamination at any wavelength (determined 
visually) from other, nearby H II regions were omitted from the 
sample. 

To correct aperture flux measurements for 'sky' contamina- 
tion, local sky levels are estimated using a centered annulus 
with an inner radius of 20' and an outer radius of 30'. The sky 
annulus is split into octants, and the median level within each 
octant is determined. In constructing the mean and RMS of 
the local sky, the highest three octants are discarded, to ensure 
against contamination of the sky level by nearby emission. In 
(rare) cases where there is no nearby emission, this approach 
results in a slight systematic underestimate of at most la in the 
sky level, or a corresponding overestimate of at most a few per- 
cent in the H II region flux. In most cases, however, this method 
largely negates the effects of emission in the sky annulus and 
reproduces manually estimated sky levels to within la with no 
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FIG. 3. — The morphology of the LMC at Ha, 1900A, R band, and in the FIR 60fim band. The 52 H II regions in our sample are circled. 



hint of systematic offset. 

Aperture corrections are typically not required for most of 
the data. However, because the FWHM of a point source in 
the 8.55 GHz radio map is 210", a modest aperture correction 
is applied. We derive this aperture correction empirically by 
smoothing the observed image with a 210" FWHM Gaussian 
beam and comparing the aperture flux in the convolved and 
original images. The ratio of these fluxes is adopted as the aper- 
ture correction. This procedure is tested by re-convolving the 
convolved image (i.e. how well can one recover the flux in the 
original image by aperture correcting the flux measured from 
the convolved image). This procedure appears to be accurate to 
~5%. Typical aperture corrections are ~ 15-25%, depending 
on aperture size and source morphology. 

Color correction in the FUV and FIR is required (the Ha 
and radio do not require color correction because of the nar- 
row bandpasses). Fluxes are corrected to the nominal effective 
wavelength of the filters. FUV color corrections, in terms of 
the observed 1500 A to 1900 A flux ratio, are calculated using 
a sequence of IUE spectra that are artificially reddened using 
a SMC extinction curve (used because it lacks a 2200A bump, 
leading to a more stable color correction). Using different IUE 
spectra and/or a LMC extinction curve to artificially redden 
the spectra leads to differences in color correction smaller than 
~3%. FIR color corrections are estimated from the observed 
60^m to lOO^m flux ratio, using the multiplicative corrections 
from the IRAS Explanatory Supplement (1.09 ±0.01 at 60^m 
and 1.01 ±0.01 at 100/im, for the range of observed Hll re- 
gion FIR colors). Color correction for the 8.3/im, 12^m, and 
25 [im fluxes was not attempted due to the complexity of the 
FIR spectrum at those wavelengths and negligible contribution 
of the FIR emission from those wavelengths to the integrated 
FIR emission. 

In Fig. 4, we show examples of our H II region SEDs (omit- 
ting the radio flux point because it makes a negligible contri- 
bution to the total bolometric flux of a star formation region, 
and the Ha data point because it is line emission). The error 
bars (which are in some cases smaller than the plotted symbol) 
reflect sky level uncertainties: a 5% uncertainty is added to the 
sky error in the FUV to account for the color correction (and 
adding an arbitrary modest contribution for 'flat fielding'). A 
5% uncertainty is added to the sky error in the radio to account 
for aperture correction errors. Ha errors include an arbitrary 
5% flat fielding and continuum subtraction uncertainty. Cal- 
ibration and image alignment uncertainties are not accounted 
for in the error bars. We do not correct for foreground galac- 
tic extinction in constructing the SEDs, but correct for it when 
quantifying the H II region attenuations. We present the H II re- 
gion fluxes and uncertainties within the 4.9' apertures in Table 
2. The fluxes in this Table have not been corrected for galactic 
foreground extinction. 

4. THE EFFECTS OF DUST IN LMC H II REGION SEDS 

Fig. 4 illustrates the significant effect of dust on H II region 
SEDs. As expected, increasing amounts of dust progressively 
depress the UV flux, re-processing this energy into the FIR. 
By using our understanding of the physics of nebular emis- 
sion and the spectra of different stellar populations, we can 
produce quantitative attenuation estimates. In this section, we 



apply the various attenuation estimators described previously 
and compare the results. For reference, our measurements of 
the H II region attenuation and UV spectral slope /3 within the 
4.9' apertures, corrected for foreground galactic extinction, are 
presented in Table 3. 

4. 1 . Comparing estimates of Ha attenuation 
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FIG. 5. — Solid points denote radio-derived Ha attenuations (Ana.Radio) 
against Balmer-derived attenuations (AH Q .Balmer), using the 4.9' aperture of 
Caplan & Deharveng (1985, 1986) for 50 of our 52 Hll regions. The Balmer- 
derived attenuation A Ha ,Balmer is derived using Ha/H/3 values from Caplan & 
Deharveng (1985). Open points denote the radio-derived attenuations from 
Caplan & Deharveng (1986), for the 29 Hll regions in common with our 
clipped sample of 50 H II regions. The solid line denotes the relationship ex- 
pected if A HctjBalmer equals A Ha ,Radio- 



Caplan & Deharveng (1986) found that radio-derived Ha at- 
tenuations of LMC H II regions, A Ha ,Radio, exceeded the Balmer- 
derived Ha attenuations, AH Q ,Baimer- This result was confirmed 
using independent radio data for LMC H II regions (Ye 1998) 
and for Hll regions in M51 (van der Hulst et al. 1988). This 
systematic offset was attributed to clumped dust, which makes 
the effective attenuation curve grayer (cf. panel (c) of Fig. 2). 

We re-examine this question using Ha/H/3 values from Ca- 
plan & Deharveng (1985) corrected for Galactic extinction, 
in conjunction with Aea.Radio values derived from our galactic 
extinction-corrected aperture Ha fluxes and aperture-corrected 
radio fluxes. We discard N30BC (no Balmer decrement) and 
N158A (radio detection less than 2a). Uncertainties in Aea.Radio 
come from propagating the uncertainties in the SED measure- 
ments and the adopted electron temperature. The uncertainty 
in A Ha! Baimer is calculated to be 9% by propagating a 4% un- 
certainty in Ha/H/3 (Caplan & Deharveng 1985) and the un- 
certainty in the electron temperature. We plot our values for 
AHa.Radio and A HQ ,Baimer, the line of equality for reference, and 
the original data from Caplan & Deharveng (1986) in Fig. 5. 
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FIG. 4. — SEDs in 4.9' apertures for four Hll regions, sorted by extinction (left to right; top to bottom): N51D (~0 mag at Ha), N113 N (~0.15 mag at Ha), 
N44D (~0.25 mag at Ha), and N159 (~0.6 mag at Ha). The Ha extinctions are determined using the Balmer decrement within 4.9' as measured by Caplan & 
Deharveng (1985), assuming a Milky Way extinction curve and a foreground screen dust distribution. The logarithm of the energy flux AF^ in WirT 2 is shown at 
1500A, 1900A, 8.3/tm, 12/im, 25/im, 60/im, and 100/tm. The 8.3/tm flux can be larger than the I2fim flux due to emission from the unidentified infrared bands, 
most plausibly associated with polycyclic aromatic hydrocarbons (e.g. Puget & Leger 1989). These data are not corrected for foreground galactic extinction, but are 
color-corrected (see text for details). Note that some error bars are smaller than the plotting symbols. 



For 29 H II regions in common we find a good overall agree- 
ment between our values of AHa,Radio and those from Caplan 
& Deharveng (1986), with a mean offset of 0.1 mag (ours are 
larger) and RMS of 0.3 mag. Omitting the 4 H II regions with 
discrepancies larger than 0.5 mag (three high values and one 
low value) reduces the offset to 0.06 mag and the RMS to 0.14 
mag. For reference, our median AH Ql R a dio formal uncertainty 
for these 29 Hll regions is 0.11 mag. We confirm the offset 
between A Ha ,Radio an d AHa.Baimer observed by Caplan & Dehar- 
veng (1986): we discuss the origin of this offset in §5.1. 

4.2. Comparing estimates offar-UV attenuation 

We now assess the consistency of indicators of the attenu- 
ation suffered by the starlight from young stellar populations. 
We compare reddening estimates, and inter-compare these with 
estimates of the attenuation in the UV derived using UV/FIR 
energy balance. When comparing these attenuation indicators, 
we also present the naive expectation from a screen model with 



LMC and SMC Bar-type dust. The intention of showing these 
models is to enable one to use existing intuition to help un- 
derstand the results. These models do not represent the true 
situation. In particular, it would be incorrect to differentiate 
between possible extinction curves on the basis of these data 
because radiative transfer effects complicate the interpretation 
of these plots. 

4.2.1. Estimating the attenuations 

In this section, we discuss the construction of FUV attenua- 
tion estimates using the FU V/FIR ratio and outline the methods 
used to estimate the stellar E(B-V) and the UV spectral slope 

13. 

As discussed in §2, it is possible to estimate FUV attenuation 
using a FUV/FIR flux ratio. This can be done using a simple 
energy balance. However, we use the flux ratio method of Gor- 
don et al. (2000) to estimate the attenuation, to allow us to fully 
take into account uncertainties due to dust geometry, type and 
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stellar populations. The 1500A attenuation derived using this 
method, as described in §2, is approximated to within 0.03 mag 
by eqn. 2. 

To apply this method we construct an integrated FIR 8- 
lOOO^m flux following Gordon et al. (2000). The 8.3/xm, 
12/im, 25/zm, 60^m and lOO^m fluxes were numerically inte- 
grated to provide an estimate of the 8 - 120yum flux. The 60^m 
and 100^m fluxes were then used to define a dust temperature 
for a (3 = 1 dust emissivity model to extrapolate to lOOO^m. 
This extrapolation, which adds typically ^25% to the flux, was 
tested against raw 178" square aperture fluxes from the ISO 
archive. To within a constant scaling that accounts for aperture 
correction and pointing offsets, the fluxes and model extrap- 
olations are within 10% of each other at 120/im, 150/im, and 
200/im. We assume no Lyman continuum ionizing radiation ra- 
diation is absorbed by dust internal to the H II region: assuming 
that 50% of the ionizing radiation is absorbed by dust produces 
around a -0.15cr systematic offset in the estimated attenuation. 

Alternatively, a more 'standard' manner to compute the at- 
tenuation of starlight from the central young cluster utilizes the 
measurement of reddening based on the optical stellar color ex- 
cess E(B-V). For 45 particular LMC Hll regions, Caplan & 
Deharveng (1986) provide E(B-V) values, which we correct 
for galactic foreground reddening. However, these values are 
average values for a number of stars near the H II region. Be- 
cause some of these stars may be well in the foreground or 
background relative to the H II region, the average reddening 
does not necessarily reflect the attenuation toward the central 
star cluster in the H II region. 

Finally, one can estimate the attenuation of the central star 
cluster using the UV spectral slope (3, as estimated using the 
1500 A - 1900 A color and as defined in eqn. 1. Uncertainties 
in f3 are derived from the random errors in the FUV flux, and 
typically become larger for more highly-reddened H II regions 
(as they are typically more attenuated and are therefore fainter). 
The (3 values of 4 UV-bright star formation regions were tested 
against 10"x20" International Ultraviolet Explorer (IUE) UV 
spectral slopes. We find that 3/4 of the IUE (3 values agreed 
with the Rocket UV-derived [3 values to better than 0.5 in UV 
spectral slope. This is a remarkable agreement considering the 
huge differences in aperture size (10"x20" compared to a 4.9' 
diameter aperture). It is worth briefly noting that for these four 
H II regions with IUE spectra there is no evidence for a system- 
atic difference between values of (3 derived using a synthetic 
1500 A - 1900 A color and using the wavelength windows de- 
fined by Calzetti et al. (1994). 

4.2.2. Comparing the attenuations 

We show the direct comparison between UV spectral slope (3 
and A i5oo.fr in Fig. 6. Despite the significant scatter, and con- 
siderable observational errors, there is a trend between (3 and 
1500 A attenuation. Two points are apparent immediately. 

Firstly, the change in FUV spectral slope (3 for a given 1500 A 
attenuation is much larger than expected using an average LMC 
extinction curve (solid line). Furthermore, the correlation is 
much steeper than the correlation between (3 and FUV attenua- 
tion derived for starbursting galaxies (open points) by Calzetti 
et al. (1994) and Meurer, Heckman, & Calzetti (1999). The 
steepness of the relationship between (3 and Aisoo.fr is likely 
due to radiative transfer effects, and does not imply that the 
dust in LMC H II regions is SMC-type: this point is discussed 
in detail in §5.2. The existence and steepness of the correlation 
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Fig. 6. — FUV spectral slope f3 against the 1500A flux ratio-derived atten- 
uation Aisoo.fr- Quantities were derived within the 4.9' aperture of Caplan 
& Deharveng (1985, 1986) for our sample of 52 Hll regions. The solid line 
denotes the expected correlation assuming a typical LMC extinction curve and 
the dashed line the SMC extinction curve. These lines, also shown in Figs. 
7-12, have been shown only to let the reader use their existing intuition to help 
understand the data; however, it is inappropriate to try differentiate between 
different extinction curves using these data because of the complicating effects 
of radiative transfer (see e.g. Figs. 2 and 12). An intrinsic fi of -2.0 has been 
assumed. Open points denote the starburst sample of Calzetti et al. (1994): 
their values of (3 were calculated in a similar way to the LMC H II regions. 
Average results from Meurer et al. (2000) for a sample of 6 ultra-luminous 
infrared galaxies are also shown, which have estimated 1500A attenuations in 
excess of 5 magnitudes. 



are robust: omitting the three points with highest Ai5oo,fr, or 
even one quarter of the data points with highest Ai 50 o.fr, does 
not significantly change either the significance or slope of the 
correlation. 

The open points were calculated using synthetic 1500 A and 
1900 A fluxes derived from IUE spectra, and FIR luminosities 
calculated as above for the H II regions using data from Calzetti 
et al. (1994, 1995). Typical errors in (3 and A 150 o,fr for the star- 
burst galaxies are ^0.1 and ~ 0.15 mag respectively, which is 
smaller than the (modest) scatter shown by the starbursts. Fi- 
nally, we also show average results from Meurer et al. (2000) 
for a sample of 6 ultra-luminous infrared galaxies that have es- 
timated 1500A attenuations in excess of 5 magnitudes. 

This comparison of empirical results demonstrates directly 
that the correlation between (3 and FUV attenuation is highly 
variable. Although H II regions are arguably the basic unit of 
star formation, the behavior of the /3-attenuation correlation 
differs strongly between H II regions and intensely star-forming 
galaxies. 

Secondly, the scatter in the Hll region /3-attenuation cor- 
relation is large. This scatter indicates that the FUV spectral 
slope is not necessarily an accurate indicator of FUV attenua- 
tion, which is not surprising because the relationship between (3 
and FUV attenuation is highly dependent on dust/star geometry 
and FUV extinction curve shape (e.g. Witt & Gordon 2000, see 
also Fig. 2). We discuss Fig. 6 in more detail in §5.2, while in- 



The Effects of Dust 



9 



vestigating possible reasons for the steep slope and scatter and 
discussing some implications of this result. 
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FIG. 7. — Optical color excess E(B-V) against the 1500A flux ratio-derived 
attenuation A1500.FR- The E(B-V) values were derived for the overall region 
around the 45 H II regions of interest. Flux ratio-derived attenuations were 
derived within a 4.9' aperture. The solid line denotes the expected correla- 
tion assuming a typical LMC extinction curve and the dashed line the SMC 
extinction curve. 
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FIG. 8. — Optical color excess E(B-V) against the UV spectral slope 
derived within 4.9' apertures. Solid and dashed lines denote the expectations 
of a screen model with LMC-type and SMC Bar-type dust respectively. 



FUV attenuation may be estimated from E(B-V), in con- 
junction with assumptions regarding the extinction curve and 



dust geometry. In Fig. 7, we show Aisoo.fr against E(B-V). 
In Fig. 8, we show the UV spectral slope (3 against E(B-V). 
Overplotted are expectations for the correlations assuming an 
average LMC extinction curve (solid line) and SMC bar ex- 
tinction curve (dashed line). Despite the large scatters in both 
correlations, the correlations are significant at greater than the 
99.9% level, with Spearman rank correlation coefficients of 
0.57 (A l50Q:FR -E(B-V)) and 0.54 (/3-E(B-V)). 

From this discussion, we conclude that the values of Ai5oo.fr, 
(3 and E(B-V) for our sample of H II regions are broadly con- 
sistent: at the very least, as one increases, the other indica- 
tors tend to increase also. The large scatter in the correlation, 
however, demonstrates that the assumptions behind each of the 
methods are continually being violated to varying degrees. 

4.3. Comparing Ha andfar-UV attenuation 




FIG. 9. — The 150()A flux ratio-derived attenuation Ai 50 o.fr against the 
Balmer-derived Ha attenuation A Ha ,Balmer. Attenuations were derived within 
the 4.9' aperture of Caplan & Deharveng (1985, 1986) for 50 of our 52 H II 
regions. Solid and dashed lines denote the expectations of a screen model with 
LMC-type and SMC Bar-type dust respectively. 

We have tested Ha and FUV attenuation estimates indepen- 
dently for internal consistency and now turn to cross-comparing 
them. Caplan & Deharveng (1986) found a highly scattered 
correlation between E(B-V) to AHa.Baimer- We explore other 
correlations between Ha and FUV attenuation indicators in 
Figs. 9 and 10, where we show the flux ratio-derived 1500A at- 
tenuation, A i5oo,fr . against the Balmer-derived Ha attenuation, 
AHa.Baimer, and the radio-derived Ha attenuation, AHa.Radio- 

In Fig. 9 we plot Ai 50 o.fr against A HQ! Baimer- Although there 
is considerable scatter, the overall correspondence between 
Aisoo.fr and A Ha ,Baimer (significance > 99.9%; Spearman rank 
correlation coefficient of 0.72) is strong. The error-weighted 
mean ratio of A1500.FR/AHa.Baimer is 3.4±0.3 with a RMS of 1.8 
for 50 H II regions. This correspondence is as expected after 
accounting for radiative transfer effects (cf. panel (c) of Fig. 2). 

In Fig. 10 we show the correlation between Aisoo.fr and 
AHa.Radio (significance > 99.9%; Spearman rank correlation co- 
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FIG. 10. — The 1500A flux ratio-derived attenuation A1500.FR against the 
radio-derived Ha attenuation Ah q .Radio- Attenuations were derived within the 
4.9' aperture of Caplan & Deharveng (1985, 1986) for 50 of our 52 H II re- 
gions. Solid and dashed lines denote the expectations of a screen model with 
LMC-type and SMC Bar-type dust respectively. 




efficient of 0.47). The error- weighted mean ratio of Ai5oo.fr/Aho: 
is 1.3 ±0.1 with a RMS of 0.7 for 50 H II regions. This relation- 
ship is more complex than the previous one between A1500.FR 
and AHa.Baimer- For example, there are a considerable number 
of H II regions with either much lower or higher A HQ) R a dio than 
would be expected on the basis of the measured A isoo.fr- This 
extra scatter will be discussed in more detail below. Despite 
some oddities, there is a basic correspondence between Ha and 
FUV attenuations in our sample of LMC H II regions, lending 
credibility to both sets of attenuation indicators. 

5. DISCUSSION 

The observed values of A HQ ,Baimer, A Ha ,R a dio, Ai 50 o,fr, P and 
E(B-V) are generally consistent, especially given the uncer- 
tainties and the range of behaviors expected from radiative 
transfer effects. However, there are two discrepancies that de- 
serve further discussion: the highly scattered correlation be- 
tween Apia.Baimer and Apia.Radio, and the deviant behavior of the 
highly scattered correlation between A1500.FR and j3. 

5.1. Are there problems with radio-based Ha attenuations? 

The systematic discrepancy between A HQ) R a dio and A Ha ,Baimer 
is well-known and generally attributed to a relatively gray effec- 
tive attenuation curve, which is due to radiative transfer effects 
(e.g. Caplan & Deharveng 1986; van der Hulst et al. 1988; Ye 
1998). We confirm the empirical result, but suggest that the in- 
terpretation should be examined in more detail. A comparison 
of Figs. 9 and 10 shows that the 1500A attenuation correlates 
with considerable scatter with Ana.Baimer, but correlates with 
a much larger scatter with AH Q .Radk>- This behavior suggests 
a limitation in the values of Ana.Radio, rather than a limitation 
with AHa.Baimer- This runs contrary to the typical argument that 
Apia.Radio is robust because it is oblivious to radiative transfer 
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FIG. 1 1 . — Distribution of the number of ionizing H II region photons, per 
unit 1500A flux. The number of Lyman continuum photons required to pro- 
duce the observed Ha flux (assuming Case B recombination) has been cor- 
rected assuming radio-derived Ha attenuations (thick solid line) and Balmer- 
derived Ha attenuations (thick dotted line). The 1500A flux has been corrected 
using the flux ratio method (the flux ratio method was earlier demonstrated 
to give attenuations which were not inconsistent with other attenuation indi- 
cators). The predicted ratio nLyc//isoo is shown for a single burst PEGASE 
Radio model with a range of ages in Myr along the arrow at the bottom. This model 
assumes 1/3 solar metallicity : models with metallicities within a factor of three 
of this value give maximum ratios less than 0.05 dex different from the 1/3 so- 
lar model. We also show the effect of slope changes of ±0.5 from the Salpeter 
value for the zero age point (error bar). Also shown is a model curve for ~ 10 5 
simulated H II regions with ages between and and 6 Myr, and total stellar 
masses between 5OOM and 8OOOM . Scatter from small number statistics 
at the upper end of the IMF and observational errors was accounted for (the 
curve for radio-derived errors is shown; the curve for Balmer-derived errors is 
essentially identical). Fluxes were derived within a 4.9' aperture. 



effects. 

We are in a position to explore the relationship between 
AHa.Baimer and AH Q ,R a dio further than existing studies because of 
our expanded wavelength coverage. Previous studies often in- 
volved only the Balmer-derived and radio-derived attenuations 
(e.g. Caplan & Deharveng 1986; van der Hulst et al. 1988; Ye 
1998), leaving the situation ambiguous. Other data (such as the 
UV and FIR fluxes) can be used, in principle, to shed new light 
on the problem. 

One sensitive diagnostic is the Ha-to-1500A flux ratio. This 
ratio is a sensitive indicator of the temperature of the illumi- 
nating stellar population. Assuming Case B recombination, 
we convert this ratio into «Lyc//i500, assuming 0.45 Ha pho- 
tons per Lyman continuum photon (Hummer & Storey 1987): 
we plot this in Fig. 11. Before converting the Ha luminosi- 
ties into the number of Lyman continuum photons, we correct 
for Ha attenuation using either radio-derived (thick solid line) 
or Balmer-derived (thick dotted line) attenuations. The frac- 
tion of Lyman continuum photons absorbed by dust before be- 
ing down-converted to Ha and radio is not accounted for. The 
loss of Lyman continuum photons to dust prior to their con- 
version to Balmer photons would simply increase the estimated 
"Lyc/ /1500: for reference, DeGioia-Eastwood (1992) estimates 
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the fraction of Lyman continuum photons lost to dust at be- 
tween 0.2 and 0.5 for six LMC H II regions (corresponding to 
an underestimate of «Lyc//i500 of between 0.1 and 0.3 dex). 
The 1500 A flux is corrected for attenuation using the flux ratio- 
derived attenuations (as is apparent from Fig. 6, attenuations 
estimated from (3 and the calibration of Calzetti et al. (1994) 
would be implausibly large and are therefore not used here). 

The radio-derived attenuations lead to larger «Lyc//i500 than 
the B aimer-derived attenuations. To be able to interpret the dif- 
ferences between the two sets of values, we show the expected 
range of «Lyc//i500 from the PEGASE stellar population synthe- 
sis model (arrow; the numbers denote H II region age in Myr). 
The illustrated model has a metallicity of 1/3 solar. Changing 
the metallicity by a factor of three in either direction, or adopt- 
ing the STARBURST99 model of Leitherer et al. (1999) results 
in systematic offsets of less than 0.05 dex. Changing the slope 
of the initial mass function from the Salpeter logarithmic slope, 
x = -1.35, by ±0.5 results in changes of ^0.1 dex in «Lyc//i500 
at zero age (the error bar in Fig. 1 1). Comparison of the model 
predictions of «Lyc/ /i500 with the observations suggests that the 
radio-derived Aea.Radio may be overcorrecting the Ha luminos- 
ity, while the Balmer-derived A Ha ,Baimer may be more appropri- 
ate. 

To test the significance of this suggestion, we construct a 
simple model. The evolution of stellar populations in narrow 
mass ranges was calculated using the PEGASE stellar popu- 
lation model. Around 10 randomly-generated stellar popu- 
lations with ages between and 6 Myr and total masses be- 
tween 500M Q and 8000M Q were extracted from a Salpeter 
IMF. Model values of «Lyc/ /1500 were convolved with the obser- 
vational errors. The model curve is shown in Fig. 1 1 . There is a 
significant tail of model log 10 (n Lyc / /1500) values above the zero- 
age value of 13.55, which are attributable in part to the effects 
of observational error, and in part to small number statistics in 
populating the IMF of H II regions with massive stars. We find 
that the observed distribution corrected with radio-derived Ha 
attenuations has a number of rather large values of «Lyc/ /1500, 
but it is not inconsistent with model distributions. The distri- 
bution corrected with Balmer-derived Ha attenuations is also 
consistent with model distributions. 

So, while there is a definite discrepancy between observed 
values of AH Q ,Baimer and A H a,Radio> it i s not possible at this stage 
to identify unambiguously the limitations of each attenuation 
indicator. Folding in FUV and FIR data can, in principle, give 
further insight, and, indeed, there is a hint of limitations in 
the radio-derived Ha attenuations. However, even given this 
extra data, it is impossible to rule out that the distribution of 
"Lye/ /1500, derived using radio-based Ha attenuations, is con- 
sistent with plausible stellar populations. More detailed inves- 
tigations, involving spatially-resolved maps of H II regions at 
Ha, H/3, and a number of radio frequencies, are required to 
fully address this question. 

5.2. The relationship between (3 and A1500.FR 

In §4.2, we discussed the highly scattered and rather steep 
correlation between the UV spectral slope (3 and the 1500 A at- 
tenuation A i5oo : fr- Both the slope and scatter of the correlation 
run counter to (at least) observational expectation: the correla- 
tion between (3 and A1500.FR for starbursting galaxies is tighter 
and shows a shallower slope. We discuss possible interpreta- 
tions of the steepness and scatter in this correlation, and some 
implications of this result for studies of the UV properties of 



high-redshift galaxies. 

We first address the possibility that the steepness and scatter 
of this correlation reflect the data quality (or, for example, varia- 
tions in Galactic foreground extinction), rather than the intrinsic 
properties of the Hll regions. The quality of the data was dis- 
cussed extensively in §§3.1 and 4.2. Briefly, the zeropoint un- 
certainties in the UV and FIR data used to construct the Aisoo.fr 
estimates are accurate to better than 10%. We have taken care to 
quantify and estimate appropriate sky level errors (§3.3), which 
are propagated into the final attenuation and UV spectral slope 
estimates. Furthermore, although the UV spectral slopes are 
difficult to check against external comparison data, we found 
an agreement of better than 0.5 in slope in 3/4 test cases (com- 
paring values derived from 10" x 20" apertures with our values 
derived using a 4.9' aperture). We conclude that our error bars 
are fair, and that data quality issues are not likely to cause the 
large scatter and steep slope we see in Fig. 6 (re-plotted as Fig. 
12). Variations in Galactic foreground extinction correspond- 
ing to AE(B-V) of ±0.06 mag produce a negligible spectral 
slope change and variations in Aisoo.fr of ±0.5 mag. This kind 
of variation is possible: Oestreicher, Gochermann, & Schmidt- 
Kaler ( 1 995) discussed the evidence for variable foreground ex- 
tinction in the direction of the LMC, but with insufficient spatial 
resolution for our purposes. Given this, we conclude that vari- 
able foreground extinction can contribute significantly to the 
scatter, but cannot account for it in its entirety. 




-4 Li , 1 , 1 , 1 , 

12 3 4 

^1 500.FR 

Fig. 12. — FUV spectral slope fi against the 1500A flux ratio-derived at- 
tenuation A i5oo.FR- Symbols are as in Fig. 6. The solid line without symbols 
denotes the expected correlation assuming a typical LMC extinction curve. 
The dashed line (which is hidden under the solid line with triangles) denotes 
the SMC extinction curve. An intrinsic f3 of —2.0 has been assumed. These 
lines have been shown only to let the reader use their existing intuition to 
help understand the data; however, it is inappropriate to try differentiate be- 
tween different extinction curves using these data because of the complicating 
effects of radiative transfer. To illustrate this point, radiative transfer model 
lines for LMC2-type dust (for a supershell near 30 Dor; triangles) and SMC 
bar-type dust (diamonds) with homogeneous (solid lines with symbols) and 
clumpy (dotted lines with symbols) geometries are overplotted. 

We turn to intrinsic properties of the H II regions to explain 
the observation. No more than 0.05 dex of scatter in (3 and 
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0.15 mag of scatter in Aisoo.fr can be generated by appealing to 
small number statistics at the upper end of the H II region IMF: 
this scatter is much smaller than the observed scatter. A promis- 
ing possibility is the combined influence of differing extinction 
curves and dust/star geometry: in §2 we demonstrated that these 
can strongly affect the relationship between Ai5oo,fr and (3. We 
demonstrate this more quantitatively in Fig. 12, where we replot 
Fig. 6, now including models from the DIRTY radiative transfer 
model. These models, in the spirit of the earlier inclusion of 
LMC and SMC screen models in earlier plots, are not meant 
to suggest that we are trying to differentiate between different 
geometries or extinction curves: rather, we are attempting to 
guide our intuition as to what effects can be plausibly caused 
by extinction curve and geometry variations. 

Variations in Hll region extinction curves and/or star/dust 
geometries can produce a steep slope and a large scatter in the 
/3-Ai5oo,fr plane. Slopes in this plot steeper than those pro- 
duced by a simple screen are possible because of scattering 
(compare the SMC Bar-type extinction curve screen and shell 
model). However, deviations from a perfect shell introduce sig- 
nificant scatter and flatten the relation somewhat (the clumpy 
shell models have shallower slopes than the homogeneous shell 
models). Lastly, it may be possible to produce redder spectral 
slopes and very low Aisoo.fr values with either Galactic fore- 
ground extinction variations or small amounts of dust directly 
along the line of sight, which produce significant reddening but 
relatively little FIR luminosity. 

We conclude that modest (and expected) variations in ge- 
ometry can produce the steep slope of the /?-Ai5oo,fr correla- 
tion and at least some of the scatter. This result carries with it 
serious implications for the use of (3 as an accurate indicator 
of FUV attenuation. Neither ultra-luminous infrared galaxies 
nor our sample of LMC H II regions follow the relatively tight 
correlation between (3 and A^oo^r established for starbursting 
galaxies. However, this correlation is utilized extensively for 
the estimation of the SFRs of high-redshift galaxies (see e.g. 
Adelberger & Steidel 2000, for extensive testing of this cor- 
relation, and references therein). Although the high-redshift 
star-forming galaxies are likely to be more similar to the low- 
redshift starburst sample of Calzetti et al. (1994) than to our 
sample of LMC H II regions, it is nevertheless a concern that 
the basic units of star formation (Hll regions) have a drasti- 
cally different behavior from star-bursting galaxies. Clearly, 
much caution is warranted before estimating FUV attenuations 
based on UV spectral slopes. 

6. SUMMARY 

We have used FUV, Ha, FIR, and 8.55 GHz radio images of 
the LMC to investigate the SEDs and attenuations of H II re- 
gions drawn from the catalogs of Caplan & Deharveng (1985, 
1986). We briefly discussed the relative morphologies of the 
LMC in the FUV, optical continuum, Ha, and FIR. In particu- 
lar, the FUV and Ha distributions within the LMC are different, 
implying that our findings regarding the effects of dust in the 
UV in LMC H II regions may not apply to the LMC as a whole. 

We constructed SEDs for the sample of H II regions. The 
effects of attenuation on the SEDs are clear: H II regions with 
larger Balmer decrements show a larger FIR to UV luminos- 
ity ratio, and typically redder UV spectral slopes (Fig. 4). To 
quantify these trends, we constructed attenuation indicators of 
both the ionized H II region gas (the Balmer decrement and ra- 
dio to Ha luminosity ratio) and the stars (the UV spectral slope 
j3, optical stellar reddening, and the balance of FUV and FIR 



luminosity). Our principal findings are as follows. 

• Estimates of Ha and FUV attenuation all correlate, although 
often with substantial scatter (Figs. 5-10). Some of this scatter 
is due to errors in the (often quite modest) attenuations. How- 
ever, much of this scatter is naturally interpreted in terms of 
radiative transfer in non-trivial dust/star geometries. 

• The modest B aimer-derived Ha attenuations are systemati- 
cally lower than radio-derived Ha attenuations. Previous stud- 
ies have attributed this systematic offset to variations in the 
extinction within the beam. In principle, more insight can be 
gained by incorporating FUV and FIR fluxes. However, the in- 
clusion of these data cannot, at this stage, unambiguously iden- 
tify the source of the discrepancy between B aimer-derived and 
radio-derived Ha attenuations. More work, involving spatially- 
resolved maps of Ha, H/3, and radio emission at a range of 
frequencies, is required to properly investigate this issue. 

• 1500 A attenuation and UV spectral slope (3 do correlate with 
each other (see Fig. 6), although with a large scatter and a sub- 
stantially steeper slope than was expected on the basis of the 
LMC extinction curve or the empirically-derived correlation 
between A 1500 and (3 for starbursting galaxies (Calzetti et al. 
1994). Ultra-luminous infrared galaxies populate yet another 
part of A1500//3 space (Meurer et al. 2000). While some scatter 
may be due to variations in foreground galactic extinction, both 
the large scatter and variation in slope between a variety of star- 
forming systems demonstrates empirically that (3 is not neces- 
sarily an accurate indicator of FUV attenuation. This point has 
been argued extensively from a theoretical perspective (Fig. 2; 
Witt & Gordon 2000). This conclusion carries with it serious 
implications for the many studies of the rest frame UV prop- 
erties of high-redshift galaxies that utilize the local correlation 
between (3 and A1500.FR for starbursting galaxies to (substan- 
tially) correct the observed UV fluxes for extinction. If, as is 
suggested by our work and the work of Meurer et al. (2000), 
the correlation between (3 andAi 5 oo,FR can have a variable slope 
and large scatter, serious thought must be devoted to the suit- 
ability of using (3 to correct UV fluxes for extinction. 
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Table 2 

h ii region fluxes through 4.9' apertures 
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71.8 ± 5.3 


0.052 ± 0.024 


N63A 


0.455 ± 0.023 


0.551 ± 0.028 


1.02 ±0.06 


1.57 ± 0.13 


1.54 ± 0.07 


6.4 ±0.1 


58.8 ± 1.2 


91.2 ± 3.4 


0.452 ± 0.024 


N59A 


0.283 ± 0.017 


0.359 ± 0.025 


5.44 ± 0.27 


7.66 ±0.14 


9.34 ±0.14 


45.9 ± 0.2 


327.7 ± 2.8 


443.1 ± 5.7 


1.376 ± 0.071 


N154A 


0.344 ± 0.023 


0.454 ± 0.037 


2.99 ± 0.24 


6.35 ± 0.42 


7.27 ± 0.41 


23.3 ± 0.3 


257.9 ± 14.6 


365.9 ± 29.8 


0.521 ± 0.045 


N64AB 


0.205 ± 0.019 


0.248 ± 0.023 


1.07 ± 0.06 


2.25 ± 0.12 


2.16 ± 0.11 


4.6 ± 0.1 


55.9 ± 2.2 


108.6 ±6.2 


0.194 ± 0.014 


N158C 


0.913 ± 0.062 


0.840 ± 0.095 


6.04 ± 0.43 


10.56 ± 0.57 


13.40 ± 0.77 


72.0 ± 1.6 


538.8 ± 33.6 


665.2 ± 54.8 


1.138 ±0.064 


N160AD 


0.262 ± 0.027 


0.435 ± 0.055 


7.36 ± 0.39 


21.07 ± 0.72 


26.03 ± 0.94 


120.3 ± 2.3 


815.5 ± 34.7 


970.3 ± 62.8 


2.418 ±0.123 


N158 


0.341 ± 0.019 


0.524 ± 0.036 


2.43 ± 0.47 


2.44 ± 0.65 


3.18 ±0.72 


15.3 ± 3.5 


140.2 ± 45.1 


151.2 ±61.8 


0.435 ±0.104 


N159 


0.099 ± 0.023 


0.271 ± 0.041 


6.21 ± 0.32 


21.44 ± 0.72 


30.73 ± 0.95 


170.0 ± 1.7 


1153.3 ±24.8 


1253.1 ±53.1 


2.808 ± 0.142 


N158A 


0.363 ± 0.026 


0.493 ± 0.040 


1.33 ±0.65 


0.16 ±0.83 


0.61 ± 1.11 


1.8 ±4.0 


10.6 ± 46.7 


1.7 ± 72.5 


0.218 ±0.137 


N160BCE 


0.307 ± 0.021 


0.439 ± 0.034 


4.51 ±0.24 


6.79 ± 0.64 


8.14 ±0.87 


36.9 ± 1.8 


369.6 ± 26.9 


526.6 ± 60.4 


1.026 ±0.051 


N175 


0.020 ± 0.003 


0.023 ± 0.003 


0.47 ± 0.03 


0.26 ± 0.75 


0.52 ± 0.98 


1.1 ±0.9 


19.2 ± 15.6 


32.8 ±46.1 


0.123 ±0.017 


N214C 


0.077 ± 0.005 


0.145 ±0.011 


2.28 ±0.11 


1.13 ± 0.17 


1.98 ±0.12 


8.2 ±0.1 


80.2 ± 0.9 


130.4 ± 3.7 


0.348 ± 0.031 


NGC 2100 


0.338 ± 0.024 


0.583 ± 0.042 


0.43 ± 0.23 


1.31 ± 0.46 


1.46 ± 0.47 


2.0 ± 1.5 


25.7 ± 23.2 


46.0 ± 46.7 


0.093 ± 0.042 


N164 


0.093 ± 0.019 


0.182 ±0.037 


2.34 ± 0.21 


2.62 ± 0.58 


3.25 ± 0.54 


14.3 ± 1.1 


163.4 ± 15.3 


210.7 ± 35.7 


0.492 ± 0.061 


N165 


0.043 ± 0.01 1 


0.083 ± 0.017 


0.54 ± 0.08 


1.05 ± 0.26 


1.12 ±0.22 


2.8 ± 0.4 


64.8 ±3.1 


95.5 ± 10.5 


0.103 ± 0.037 


N163 


0.019 ± 0.007 


0.056 ± 0.010 


1.27 ±0.14 


3.16 ±0.29 


3.82 ± 0.30 


8.6 ± 0.6 


93.5 ± 10.4 


160.9 ± 23.7 


0.227 ± 0.046 


N180AB 


0.357 ± 0.018 


0.515 ± 0.026 


3.44 ±0.17 


2.90 ±0.11 


3.19 ± 0.22 


13.3 ± 0.2 


120.6 ± 2.6 


199.2 ±11.3 


0.577 ± 0.032 
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Table 3 

h ii region attenuation estimates and uv spectral slopes through 4.9' apertures 



Name 



/in 



(mag) 



(mag) 



0.31 ± 0.09 
0.12 ± 0.09 
0.47 ± 0.09 
0.29 ± 0.09 
0.39 ± 0.09 
0.46 ± 0.09 
0.07 ± 0.09 
0.20 ± 0.09 
0.24 ± 0.09 
0.10 ±0.09 
0.20 ± 0.09 
0.15 ±0.09 
0.12 ±0.09 
0.15 ±0.09 
0.06 ± 0.09 
0.14 ±0.09 

-0.03 ± 0.09 
0.05 ± 0.09 
0.18 ±0.09 
0.05 ± 0.09 
0.25 ± 0.09 
0.25 ± 0.09 
0.46 ± 0.09 
-0.07 ± 0.09 
0.17 ± 0.09 
0.05 ± 0.09 
0.27 ± 0.09 
0.02 ± 0.09 
0.10 ±0.09 
0.19 ±0.09 
0.11 ±0.09 
0.10 ±0.09 
0.23 ± 0.09 
-0.06 ± 0.09 
0.08 ± 0.09 
0.50 ± 0.09 
0.19 ±0.09 
0.11 ±0.09 
0.37 ± 0.09 
0.44 ± 0.09 
0.19 ±0.09 
0.65 ± 0.09 
0.16 ± 0.09 
0.27 ± 0.09 
0.31 ± 0.09 
0.15 ±0.09 
0.08 ± 0.09 
0.50 ± 0.09 
0.25 ± 0.09 
0.53 ± 0.09 
0.16 ±0.09 



0.62 ±0.11 

1.48 ± 0.11 
1.26 ±0.11 
0.64 ±0.11 
0.57 ± 0.30 
0.74 ±0.11 
0.54 ±0.14 
0.88 ±0.10 
0.53 ±0.11 
0.68 ±0.11 
0.75 ±0.15 
0.75 ±0.11 
0.67 ±0.11 
0.86 ±0.11 
0.78 ±0.11 
0.29 ±0.13 
-0.18 ± 0.15 
0.32 ±0.11 
0.76 ±0.11 
0.80 ±0.11 
0.65 ±0.15 
0.78 ±0.11 
0.17 ±0.17 
1.30 ±0.12 
0.26 ±0.13 
0.53 ±0.14 
0.60 ± 0.27 
0.28 ±0.11 
0.19 ±0.14 
0.34 ±0.14 
0.54 ±0.13 
0.38 ±0.11 
0.02 ±0.13 
0.34 ±0.11 
-0.41 ± 0.51 

1.49 ±0.11 
0.89 ±0.11 
0.49 ±0.15 
0.53 ±0.12 
0.57 ±0.12 
1.17 ±0.11 
0.51 ± 0.34 
1.52 ±0.11 
0.42 ± 0.87 
0.77 ±0.11 
0.94 ±0.18 
0.34 ±0.13 
0.73 ± 0.76 
0.69 ±0.18 
0.58 ± 0.42 
0.51 ± 0.26 
0.44 ±0.11 



(mag) 



1.17 ± 0.21 
1.55 ± 0.32 

1.13 ±0.17 
0.69 ±0.16 
1.08 ± 0.25 
0.96 ±0.17 
0.09 ± 0.05 
0.49 ±0.13 
0.31 ±0.09 
0.64 ±0.15 
2.01 ± 0.41 
0.08 ± 0.05 
0.06 ± 0.05 
0.84 ±0.17 
0.56 ±0.14 
0.20 ± 0.07 
0.63 ± 0.15 
0.10 ± 0.05 
0.46 ±0.12 
0.37 ±0.11 
0.44 ±0.12 

1.23 ±0.17 
0.94 ±0.18 
0.86 ± 0.17 
0.02 ± 0.04 
0.32 ±0.10 
0.17 ±0.07 
0.13 ±0.06 
0.13 ±0.06 
0.13 ±0.06 
0.44 ± 0.12 
0.30 ± 0.09 
0.22 ± 0.08 
0.19 ± 0.07 
0.23 ± 0.09 
0.18 ± 0.07 
1.15 ±0.17 
0.85 ±0.17 
0.43 ±0.13 
0.71 ±0.16 
1.90 ± 0.26 
0.48 ±0.13 
3.11 ±0.42 
0.01 ± 0.04 

1.14 ±0.17 
1.04 ± 0.20 
1.10 ±0.17 
0.11 ±0.06 
1.34 ± 0.26 

1.24 ± 0.29 
2.43 ± 0.52 
0.48 ±0.12 



Ay 

(mag) 



0.72 ± 0.20 
1.28 ± 0.30 
0.52 ±0.16 
0.39 ±0.13 
0.57 ± 0.22 
0.58 ±0.17 
0.07 ± 0.05 
0.37 ±0.12 
0.21 ± 0.09 
0.42 ±0.13 
1.06 ± 0.20 
0.05 ± 0.05 
0.06 ± 0.05 
0.64 ±0.17 
0.45 ±0.14 
0.14 ±0.07 
0.34 ±0.12 
0.13 ± 0.07 
0.38 ± 0.13 
0.39 ±0.13 
0.30 ±0.11 
0.67 ±0.18 
0.62 ±0.17 
0.58 ± 0.16 
0.04 ± 0.04 
0.18 ±0.08 
0.14 ±0.08 
0.20 ± 0.08 
0.14 ± 0.07 
0.14 ±0.07 
0.30 ±0.11 
0.23 ± 0.09 
0.20 ± 0.08 
0.21 ± 0.09 
0.13 ±0.08 
0.14 ±0.07 
0.94 ± 0.20 
0.64 ±0.18 
0.33 ±0.12 
0.70 ± 0.20 
1.41 ± 0.24 
0.30 ±0.11 
2.10 ± 0.30 
0.01 ± 0.04 
0.86 ± 0.20 
0.87 ± 0.22 
0.66 ±0.18 
0.06 ± 0.05 
0.86 ± 0.25 
0.77 ± 0.24 
1.40 ± 0.28 
0.32 ±0.11 



0.48 ± 0.72 
-0.92 ± 1.28 
1.85 ± 0.36 
0.24 ± 0.41 
1.07 ± 1.25 
0.25 ± 0.35 
-1.11 ± 0.28 
-1.17 ±0.28 
-0.63 ± 0.28 
-0.46 ± 0.30 
2.39 ± 1.15 
-0.27 ± 0.28 
-1.44 ±0.28 
-0.95 ± 0.28 
-1.27 ±0.31 
-0.72 ± 0.28 
0.42 ± 0.28 
-2.83 ± 0.29 
-1.48 ±0.46 
-2.40 ± 0.28 
-0.67 ± 0.28 
1.21 ± 0.29 
-0.19 ±0.40 
-0.35 ±0.31 
-2.89 ± 0.28 
-0.00 ± 0.36 
-1.43 ±0.70 
-3.52 ± 0.29 
-2.22 ± 0.29 
-2.40 ± 0.32 
-0.75 ± 0.41 
-1.09 ±0.30 
-1.75 ±0.28 
-2.54 ± 0.30 
-0.08 ± 0.95 
-1.23 ±0.28 
-1.05 ±0.36 
-0.89 ± 0.41 
-1.24 ±0.51 
-2.30 ±0.51 
-0.00 ± 0.63 
-0.30 ± 0.35 
1.92 ± 1.06 
-0.78 ± 0.42 
-0.59 ± 0.40 
-1.31 ±0.69 
0.50 ± 0.39 
0.14 ± 0.39 
0.61 ± 1.11 
0.60 ± 1.32 
2.24 ± 1.62 
-0.55 ± 0.28 
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